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Abstract: According to general relativity, the frequency of electro- 
magnetic radiation is altered by the expansion of the universe. This 
effect — commonly referred to as the cosmological redshift — is of utmost 
importance for observations in cosmology. Here we show that this redshift 
can be reproduced on a much smaller scale using an optical analogue inside 
a dielectric metamaterial with time-dependent material parameters. To this 
aim, we apply the framework of transformation optics to the Robertson- 
Walker metric. We demonstrate theoretically how perfect redshifting or 
blueshifting of an electromagnetic wave can be achieved without the 
creation of sidebands with a device of finite length. 
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1. Introduction 

Recently, there has been considerable interest in developing optical analogues of black holes in 
table-top devices HI |2l S SJ . The electromagnetic parameters of these devices are designed in 
such a way as to mimic effects arising in the presence of black holes. Using this approach, an 
artificial event horizon has been implemented by means of a moving medium I!], an inhomoge- 
neous refractive index profile was proposed to mimic celestial mechanics 121, and a broadband 
omnidirectional absorber has been developed mi?). Although the equivalence between gravita- 
tional effects and optical media was suggested by several authors ||5]|6l|71, it has recently been 
conceived from a different point of view JS) and has become a tool to design optical devices 
within the framework of transformation optics ll9l [T0lfmfT2l . 

Transformation optics provides a new way of looking at the interaction between light and 
matter. This recently developed method in electromagnetism is based on the analogy between 
the macroscopic Maxwell's equations in complex materials and the free-space Maxwell's equa- 
tions on the background of an arbitrary metric. It describes how a general topological deforma- 
tion of reality can be implemented using materials with nontrivial constitutive parameters. The 
invisibility cloak, which bends the electromagnetic reality, i.e., the electromagnetic energy flow, 
around a hole — which therefore turns invisible — , is probably the most famous example of this 
technique ||9] [lOl [T3] [T4ll . Moreover, the cloaking principle extends beyond electromagnetism 
and can be applied in acoustics ifTSll . hydrodynamics (TSll, and even quantum mechanics ifTTll . 

Up to now, most applications of transformation optics have been related to the spatial ma- 
nipulation of optical beams and the design of imaging systems lfT8l[T9ll20ll . These devices are 
generated through spatial coordinate transformations and do not consider curved geometries. 
Transformation optics, however, applies to general four-dimensional metrics, and thus extends 
beyond the scope of existing applications. In this paper, we show how a four-dimensional met- 
ric that plays a role in general relativity — the Robertson-Walker metric — can be implemented 
directly with dielectrics and we demonstrate theoretically that the cosmological redshift can be 
reproduced inside a dielectric structure of finite size with time-varying electromagnetic proper- 
ties. 




Fig. 1. A graphical representation of a traveling wave solution inside a medium that is 
the electromagnetic analogue of the Robertson- Walker metric, (a) The spatial variation at 
constant time t = I and (b) the temporal variation of this solution for fixed location z = 1, 
when we modulate the permittivity and permeability with a{t) = l+t. 



2. Electromagnetic analogue of the cosmological redshift 

The cosmological redshift of an electromagnetic wave traveling from one place to another is an 
effect that originates from the expansion of the universe. Such an expanding background can 
be modeled by the Robertson- Walker metric. 



ds^ = -c^-dt^+a-{t) 



'''' '-da' 



(1) 



where a{t) is a real, dimensionless scale factor, related to the well-known Hubble parameter by 
H = d/a, and K represents the spatial curvature of the metric. An electromagnetic wave with 
frequency OJem emitted at instance t = fem will be observed at a different frequency Oobs at time 
t = fobs- These two frequencies are related by (see, e.g., Ref. lISTl ) 

In an expanding universe, the frequency of electromagnetic waves is redshifted. 

To construct an electromagnetic analogue of this effect, we proceed by calculating the ma- 
terial properties of a dielectric medium using transformation optics ||9] [H] [l2l . It is straight- 
forward to show that the Robertson- Walker metric without spatial curvature (fc = 0) can be 
translated into an isotropic, homogeneous dielectric with constitutive parameters 

ex = M"x = «(0, e-v = At-v = fl(0, £\ = l^\ = a{t). (3) 

The traveling wave solutions inside such a medium are found by solving Maxwell's equations 
in the absence of free sources and currents. In combination with Eq. (O, Maxwell's equations 
combine into a generalized wave equation, 



where c is the speed of light in vacuum. Its solutions are waves with time-dependent amplitude 
propagating along an arbitrary direction labeled by the unit vector l^: 

E{r,t)^-LG(r-K-c f (5) 

where G{x) is an arbitrary differentiable function. To gain some insight in these traveling waves, 
we consider a harmonic function G{x) = cos(jc) in Fig. [Hand we plot the temporal and spatial 



a) (11) (III) 



a(t) 




Fig. 2. A setup where we implement a finite Robertson-Walker device. At the leftmost 
boundary a monochromatic wave of frequency coq and amplitude Aq illuminates the device. 
We modulate the permittivity and permeability as e(t) = il{t) = a{t) = l+t. We show that 
the output undergoes a frequency shift according to the cosmological redshift formula. 

variation of Eq. (|5]l. We see that the wave has a well-defined wavelength, but that its frequency 
is changing in time. A straightforward calculation of the instantaneous frequency COinstlO = 
— (3r(k-r — \k\cJi^dt/a{t)) reveals that the frequencies at two different times fi and t2 satisfy 
Oinst(f2)/'Binst(fi) = <^{ti) / o{t2), in agreement with the redshift formula of Eq. 

3. A frequency converter of finite size 

As shown above, a traveling wave in an infinite Robertson-Walker medium is subject to an 
equivalent of the cosmological redshift. This does not guarantee, however, that we can use this 
medium to alter the frequency of an incident wave, since any practical application requires 
materials of finite extent. Furthermore, as was pointed out in Ref. 1221 . the boundary conditions 
play a crucial role in a transformation-optical device of finite extent. We therefore examine 
in this section whether a slab of the Robertson- Walker medium with thickness L (see Fig. |2]l 
would suit as a frequency converter. To simplify notation, we introduce the function 

m= / (6) 

Jto a{t) 

We illuminate the dielectric from the left in region (I) with a monochromatic plane wave 
of frequency coo = |ko|c, and we want to calculate the wave that is emitted in region (III). 
In general, each interface would give rise to transmitted and reflected waves, and one would 
explicitly impose continuity of the tangential components of both E and H. However, since the 
impedance rj = {/i/eY^^ is the same on both sides of each interface, we can anticipate that 
there are no reflected waves. We therefore restrict our attention to the electric field E, consider 
purely right-moving waves and impose continuity at each interface. In regions (I), (II), and 
(III), respectively, the electric field can be written as 



Ei(z,f) = Aoe'^of-^-^-^'lE, (7) 

Ett(z,0 = fcf^le, (8) 
a{t) 

Eiii(z,f) = //(z-cf)lE, (9) 
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Fig. 3. The wave fronts inside the frequency converter as a function of space and time when 
the leftmost boundary is illuminated by a monochromatic wave. In this diagram, we have 
simulated a linear evolution of the scale factor a{t) = t. 



where G{x) and H{x) are differentiable functions. The continuity of the tangential component 
of the electric field at z = yields 

Aoe-V' = £(z^. (10) 

This equation can be used to evaluate G{x) for an arbitrary value of its argument x: 

G(x)=afr'f--^))Aoe-'^-/-'(-f), (11) 



where we have introduced the inverse function /^'. In Fig. [3] we plot the evolution of the 
wavefronts inside the material as a function of space and time. We can now combine Eq. ^ 
with Eq. ( fTTT i to determine the electric field at the rightmost boundary [z = L) of the device: 

E„(L,0 = "(/''(/(;)-§)) ^^e-^o.:/-'(/(0-§)lB. (12) 
a{t) 

To retrieve a general expression for the electric field in region (III), we have to express the 
continuity of the electric field at the rightmost boundary and we ultimately find that 

Em(z,f)-Eii(^L,r-^y (13) 

The instantaneous frequency of the electromagnetic wave emitted in the vacuum region (III) 
behind the device {z — L) can be written as 

a(t) 

This equation shows that the proposed design from Fig. |2] indeed executes a frequency shift. 
Since the velocity of the wavefronts v = c/ ,JeJl ~ c/a{t), we know that L ~ Jl^ c / a{t) dt = 
cifih) — f{h)), where t\ and t2 indicate the time of incidence and departure of a wavefront. 
If we are observing at time t at the rightmost boundary, then (/(O ~L/c) corresponds 



to the time when the wavefront was at the leftmost boundary. Hence, we retrieve the cosmo- 
logical redshift formula, fi)inst(f2)/fi^nst(fi) = a{t{)/a{t2)- The reader should also note that the 
amplitude is altered with a similar factor, so that the device might also be used as an opti- 
cal amplifier It will often be beneficial to create a constant frequency shift, i.e., to generate a 
monochromatic wave in region (III); in this case, the right hand side of Eq. ( fT4b has to be con- 
stant. By renaming /(f) =x, this condition immediately implies that a{f^^{x)) ~ exp(ax)p(x), 
where a is an arbitrary (real) constant and p{x) is a periodic function with period L/c. Using 
Eq. ©J we observe that = and infer that the condition is satisfied when 

(x) = exp(ax)q'(x) — j3 with j3 constant and q{x) periodic with period L/c. By inverting this 
relation to obtain f{t) and differentiating with respect to t, we can calculate those a{t) which 
generate a monochromatic wave in region (III). For instance, the special case q{x) = 1 leads to 
/(f) = ln(f + j3)/a and thus to the linear profile a{t) = a{t + ji). The resulting frequency shift 
is then given by 

(OO 

4. Conclusion and discussion 

We note that the output signal has all its energy at the desired frequency without the creation of 
sidebands. The conversion is thus characterized by a very high efficiency. The proposed device 
achieves a frequency shift using linear materials, so that the superposition principle remains 
valid and an arbitrary wavepacket can be frequency-shifted. This frequency shift is possible 
due to the time evolution of its parameters, which renders it non-stationary ||23l . 

Equation (fTSl ) enables us to calculate the characteristics of a device that generates the desired 
frequency shift of a monochromatic wave. Let us suppose, for instance, that we want to shift an 
incident frequency of /o = 100 THz by 10 GHz in a material of thickness L = 2Ao = 6 \xm. We 
then find that the modulation rate of the refractive index must equal a w 5 GHz. Since the corre- 
sponding wavelength of this modulating signal is much larger than the dimensions of the device, 
we can generate the permittivity variation by low-frequency electro-optic modulation. The most 
difficult part is the time evolution of the permeability. This can be achieved by introducing an 
array of split-ring resonators 1241 l25l in the electro-optic material. The frequency-dependent 
permeability of an array of split rings can be modeled by 

Fco^ 

^{(o) = l + — (16) 

where F is the filling factor and Olc is the resonance frequency, which is inversely proportional 
to the square root of the material's permittivity ec inside the gap of the split rings. Operating far 
below resonance, we can approximate the permeabiUty by /x(a)) sa 1 + F(0^c^l^ec/{dw), with 
/, d, w the geometrical parameters of the split rings ll25l . By properly choosing the electro-optic 
material inside the gap of the split rings, we can thus vary the permeability and permittivity 
according to Eq. Q with electro-optics to establish a frequency shift determined by Eq. ( fTST l. 
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